function [x, q] = mptm(x0, q0, a, b, c, k)
% MPTM 模型函数
%   [x_next, q_next] = mptm_model(x, q, a, b, c, k) 根据当前状态 (x, q) 和参数
%   计算下一时刻的状态，模型为：
%
%     x_next = a*(x + b*sin(3*x) + c)^2 + k*cos(q)*x
%     q_next = q - x
%
% 输入：
%   x, q  - 当前状态变量
%   a, b, c, k - 模型参数
%
% 输出：
%   x, q - 下一时刻状态

    x = a*( x0 + b*sin(3*x0) + c )^2 + k*cos(q0)*x0;
    q = q0 - x0;
end
